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We review the attractor properties of the simplest chaotic model of inflation, in which a 
minimally coupled scalar field is endowed with a quadratic scalar potential. The equations 
of motion in a flat Friedmann- Robertson- Walker universe are written as an autonomous 
system of equations, and the solutions of physical interest appear as critical points. This 
new formalism is then applied to the study of inflation dynamics, in which we can go 
beyond the known slow-roll approximation. 
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1. Introduction 

One of the most studied issues in Cosmology is the dynamics of cosmological scalar 
fields, mostly because of their usefulness in providing models for different needed 
processes in the evolution of the universe. It is not an easy task at all to present a 
full bibliography on scalar fields, but a nice and comprehensive review was recently 
presented by Copeland, Sami and Tsujikawa in Ref |TJ 

A full plethora of methods exist in the specialized literature to study the dy- 
namics of cosmological scalar fields, but only a few can be told to be extreme useful 
and of general applicability, in a field where exact solutions are rarely found. 

One of these methods is the writing of the scalar and gravitational evolution 
equations in the form of a dynamical system. This is an idea that has largely per- 

O Q A 

vaded the specialized literature on cosmological scalar fields . The reason is that 
the theory of dynamical systems can show the existence of (fixed) stationary points 
that may represent important cosmological solutions; many times these solutions 
are attractors the system evolves to independently of the initial conditions. This is 
of wide interest in the case of inflationary and dark energy modelJ^^. However, the 
dynamical system approach is not always completely suitable for the study of scalar 
field dynamics. 

A typical example is a scalar field model endowed with a quadratic potential. 
The first attempt for this case was made by Belinsky et aJ2i (see alsd^, in which 
the dynamical variables were just the normalized values of the scalar field and grav- 
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itational variables (the scalar field value (f), its time derivative 0, and the Hubble 
parameter H). Their analysis did not reveal inflationary attractor solutions, but 
rather asymptotic inflationary behavior at times t ^ — oo, for which the corre- 
sponding fixed points were all unstable. 

In a later paper, de la Macorra and Piccinelll^ used the same variables suggested 
in Ref.[3|for an exponential potential; for the latter is possible to write the equations 
of motion as an autonomous plane system. On the contrary, the system of equations 
for a quadratic potential case cannot be written in such a simple form, and one 
has to introduce extra variables, which are not directly related with the original 
dynamical variables, to have a closed system of equations. In any case, de la Macorra 
and Piccinelli were able to obtain asymptotic solutions, and provided a general 
classification for the behavior of arbitrary scalar potentialJ^. 

Our own proposal to study the dynamics of a quadratic potential is so meho w in 
between the last two discussed above. We take the same variables used i and 
define a third new variable related to the Hubble parameter; the resulting dynamical 
system is a 3-dimensional autonomous one. However, the nature of the fixed points of 
physical interest is clearly revealed if one studies a reduced 2-dimensional dynamical 
system, and the third variable is just taken as a control parameter. 

The critical points that will be considered are not fixed in the strict sense, as 
their location on the 2-dimensional phase space depend upon the values of the 
control variable. Nevertheless, we shall show that the critical points of the reduced 
dynamical system carry useful physical information, and that their stability can be 
established by standard methods. 

A similar approach was presented in a recent paper irP, in which the dynamics 
of quintessence models is under the control of a so-called roll parameter A oc V /V. 
The location of the fixed points, which are actually those found for an exponential 
potentiaPJ, is determined by the value of A at any instant of time. Such a property 
provides a powerful method to study general quintessence models. 

Even though a qua dratic p otential is the simplest and the most well known of 
the infiationary 

niodeiPsnnr 

we will revise some inflationary solutions from the 
point of view of the results on attractor solutions obtained from the dynamical 
system approach. 

A summary of the paper is as follows. In Sec. [21 we describe the mathematical 
definitions of the new dynamical system, and make a detailed study of its critical 
points. In Sec.[3l we apply our results to inflationary models in the simplest chaotic 
scenario. Finally, conclusions are discussed in Sec. HI 



2. Mathematical background 

The model is defined by the following action in the scalar field 



(1) 
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where m indicates the scalar field mass. The evolution equations in a spatially flat 
FRW model, with a metric given by gp,,, = diag(— 1, a^(i), a^(<), a^(t)), are 

= -Sff^-TO^^, (2a) 
H = -AttG^^ , (2b) 

together with the Friedmann constraint 



H 



^ 2 

2 



In all the above equations, a dot means derivative with respect to the cosmic time 
t, a{t) is the scale factor, H is the Hubble parameter, G is Newton's gravitational 
constant, and we use units in which c — 1. 



2.1. Dynamical system 

To study the dynamics of the system ([2]), it is convenient to define new dimensionless 
variables as follows. 



VIttG ■ VIttG , m 

a; = — = — ffl , y = — = — m0, z= — . (4) 

The evolution equations ^ are then replaced by the equations 

X — —3(1 — x'^)x — zy , (5a) 
y' — 3x^y + zx , (5b) 
z' = 3x^z . (5c) 

where a prime denotes derivative with respect to the e-fold number iV = Ina, and 
we have used Eq. ([3]) in the form —H/H^ = 3x^; on the other hand, the Friedmann 
constraint now reads 



x 



2 



y' = i. (6) 



Eqs. ([5]) are written in the form of a dynamical system, which helps to clarify the 
properties of the original equations of motion. 

Unlike the very well known case of an exponential potentiaPEd, it is not 
possible to take the Friedmann parameter H out of the evolution equations, and we 
are forced to take it into account through the definition of the new variable z; not ice 
that this new variable is, by definition, a monotonic growing function (seil2Il31 for 
a similar situation in a cosmological dynamical system). 

The trajectories in the 3-dim phase space are located on the surface of an infinite 
cylinder of unitary radius, and, strictly speaking, the only critical point is the origin 
of coordinates. It is then evident that the 3-dimensional system does not provide 
enough relevant information about the dynamics of the original physical system. 

However, one is usually interested in the evolution of variables x and y, whose 
quadratic values represent the contributions of the kinetic and potential energies. 
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respectively, to the total energy contents of the universe. Because of this, the Fried- 
mann constraint involves the values of x and y only, see Eq. ([3]) , and then all im- 
portant physical features of the original system of equations ^ are more tractable 
if we restrict ourselves to the xy plane of the phase space. 

Thus, we are to study the reduced 2-dim dynamical system provided by Eqs. ((5a|) 
and (jSbp only; and variable z will be considered a control parameter all the results 
presented below will depend upon. In what follows, and for purposes of generality, 
we will study the unconstrained system (|5ap and (jSbp . and we shall impose the 
Friedmann constraint ^ whenever is necessary to clarify the nature of the critical 
points. 

2.2. Critical points 

The critical points of the 2-dimensional system correspond to those for which x' = 
y' = 0. The results can be summarized as follows. 

Null scalar field. This point corresponds to the origin of coordinates x ^ y = 0. 
This point exists for any value of z; but it should be stressed out that this point 
does not accomplish with the Friedmann constraint ([6]), i.e., it is inaccessible to the 
original physical system. 

Scalar field domination. There are four critical points with coordinates given by 



Notice that these critical points cannot be considered fixed in the strict sense, as 
their values depend upon variable z; but it is remarkable that they all lie, whenever 
their existance is allowed, on the unitary circumference Xg + yg = 1. 

For z = 0, the four critical points are located at {±1,0} and {0,±1}; but 
as z evolves to larger values, the points move one to each other in the unitary 
circumference until they merge into two points located at {±l/-\/2, =f1/-\/2} once 
z — 3/2. The critical points cease to exist for z > 3/2, see also Fig.[T] 

2.3. Stability analysis 

We now proceed to the stability analysis under the assumption that variable z is, as 
mentioned before, only a control parameter. Following standard methods, we make 
small perturbations around each one of the critical points of the form x = xg + u, 
where x = {x,y), u = [u^, Uy) are the corresponding perturbations of the dynamical 
variables, and a subscript '0' denotes the critical points. 

The dynamical system in Eqs. (|5ap and (|5b[) is of the form x' = f(x, z), and 
then the evolution equations for the perturbations are u' = A^u, where the stability 
matrix M has components 
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We only need to calculate the eigenvalues of the stability matrix by taking the values 
of xq and at each one of the critical points; the results are summarized below. 
Null scalar field. The eigenvalues of the stability matrix in this case are 

(9) 

For any value z > 0, this critical point is stable; for values < z < 3/2, the point is a 
node, whereas for values 3/2 < z the point is a converging focus. Unfortunately, we 
have already mentioned that this critical point does not accomplish the Friedmann 
constraint ©J and then it is not of physical interest. 

Scalar field domination. For the rest of the critical points ([7]), the eigenvalues 
are given by 

UJl = 6xg , tJ2 = 3 {2x1 - 1) ■ (10) 

It can be seen that all the critical points are unstable under small perturbations 
whenever they exist, as the first eigenvalue is always positive, oji > 0. However, 
the second eigenvalue is negative (positive) for Xq < 1/2 {xq > 1/2). This change 
of sign does not make any difference for the stability of the critical points, except 
in the case of trajectories subjected to the Friedmann constraint ([6]), as we are to 
show now. 

The general solution to the perturbations is of the form 

u = Ci ( """^ e^^' + C2 ( y° )e^-\ (11) 

where the terms in brackets are the normalized eigenvectors of the stability matrix, 
and the C's are arbitrary constants which are determined from initial conditions. 
The Friedmann constraint ([6]), to first order in perturbations, reads 

xou^ + youy = Cie'^'' ^ , (12) 

and then Ci = at all times. In other words, the Friedmann constraint docs not 
allow the existence of radial perturbations. 

Therefore, the nature of the critical points for any trajectory constrained to the 
unitary circumference only depends on the value of W2- In consequence, there are 
two stable (attractor) points at 




and there are two unstable points at 

Fig-fflshow some sketches of the phase portrait, as depicted in the mathematical 
package MAPLE, of our dynamical system on the xy plane for different values of 
z; the figures speak by themselves and illustrate well the description of the critical 
points made in the text above. 
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(a)z = 0.5 




(b)^ 





(c)z = 1.5 



(d)z 



Fig. 1. 2-dimensional phase portraits for the set of Eqs. llSal l and llSbl l. in which variable z is taken 
as a free parameter. The arrows show the normahzed velocity field with components v = {x',y'), 
and the unitary circumpherence represents the Friedmann constraint l(6)l. As expected, the origin 
of coordinates is a convergent focus in the four figures, and there are other four critical points 
given by Eq. (0, which are all unstable, located on the unitary circumference for < z < 3/2; the 
latter merge into two points at the special value z = 3/2. See the text for more details. 



After the attractor points has disappeared, the trajectories move in the coun- 
terclockwise direction with increasing angular frequency. To leading order for large 
values of z, one can show that variables x and y obey a harmonic oscillator equa- 
tion where variable z plays the role of the angular frequency. This fact makes the 
numerical evolution very difficult to follow for times corresponding to z > 2/3. 
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2.4. 'Fixed point' (FP) approximation 

The method that was developed in the previous sections is going to be called the 
'fixed point' (FP) approximation. The reason is that the 2-dim system (|5ap and (|5bp 
is not an autonomous one, and then the values (x, y) for which the terms on their 
r.h.s. vanish are not, in the strict sense, solutions of the 2-dim dynamical system. 

However, the 2-dim velocity field vanishes at the critical points, see Fig. [TJ 
which indicates that they should have attractor properties. In this respect, the 
standard perturbation analysis gives correct information about the stability of the 
critical points. It should be stressed out that the stability analysis is valid only on 
constant-z slices of the 3-dim phase space; in other words, we have to neglect the 
time evolution of variable z. In consequence, the FP is expected to work well only 
in the cases for which the value of z', see Eq. ([5c| . is small. 

An example is shown in Fig. [2] to illustrate that the FP approximation is very 
good at early times, and that it differs from the true evolution of the system at 
late times when z' cannot be neglected. In any case, in inflationary studies one is 
interested in the early evolution of a given model where the FP approximation is 
particularly useful. 

To finish this section, we would like to compare our method with that used in 
Ref. 121 where infiationary solutions were shown to be characteristic of the equations 
of motion Using other dynamical variables different to ours in Eqs. ([U, the 
infiationary solutions were identified with two horizontal separatrices on a 2-dim 
phase space, see Fig. 1 in Ref. [51 We show in our Fig. [21 that those inflationary sep- 
aratrices are well described by our stable critical points (|13p . From the comparison 
above, we see that the FP approximation may provide more semi-analytical results 
than other methods. 

3. Inflationary dynamics 

One of the issues cosmologists are most interested in are the infiationary solutions 
provided by scalar field models. The minimally-coupled scalar field endowed with a 
quadratic potential has a long tradi tion as an infiationary model, and its predictions 
are very well known, see 

forinstancPaiMM 

besides, it seems to be the best model 

to fit the available observational 

datiPm, see alscP. In this section, we will rewrite 
part of the inflation formalism in terms of our new variables, and discuss the role 
played by the critical points studied in the previous section. 

3.1. Attractor inflationary solutions 

The condition for an accelerated expansion is 

-=H + H^>0 ^ 3x2 = -— <1. 
a 

The description of the dynamics of a scalar field with a quadratic potential we 
provided in Sec. [2l confirms the existence of attractor points that provides an accel- 



April 17, 2009 18:13 WSPC/INSTRUCTION FILE quadratic 



8 L. Arturo Urena-Lopez and Mayra J. Reyes- Ibarra 





Fig. 2. (Top) Comparison between the numerical solutions of the full 3-dim equations of mo- 
tion JSJ (red solid lines) and the fixed point (FP) (blue dotted line) and slow-roll (SR) (green 
dashed line) approximations. The numerical solutions were given initial conditions according to 
the Friedmann constraint ^ at zi = 0.08. The trajectories for the FP approximation were ob- 
tained from the attractor solutions JTSj, whereas those corresponding to the SR approximation are 
given by Eqs. I IA.21 1 in [Appendix A| Initially, the numerical solution quickly reaches the attractor 
solution given by the FP approximation, but they differ at late times. However, the separation 
from the SR approximation is much larger; see [Appendix A| for more details. (Bottom) 2-dim 
phase space with the same variables used in Fig. 1 of Ref. 2. The red solid lines represent difi^erent 
solutions of the equations of motion JSj, whereas the green dashed lines are the trajectories ob- 
tained from the stable critical points l|13| l. Notice that the latter agree well with the inflationary 
separatrices described in Ref. 2. There is a small discrepancy close to the origin, which shows the 
break-down of the FP approximation at the late stages of inflation; see the text for more details. 
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crated expansion; but that happens only as long as the attractor points accomplish 
with the condition (fTSl) : that is, if 3a;^ < 1. 

As the position of the critical points depends upon z, we see that the attractor 
points (|13p become non-inflationary for z > Zend = V^, where the subscript end 
means the time the attractor solution leaves the accelerated regime; the correspond- 
ing values for the other variables are = 1/3 and end = 2/3. 

In fact, our calculations show that the scalar field value once the attractor is no 
longer inflationary is 

Vend mpi 
Zend ^V'"' 

where mpi = G~^^^ is the Planck mass. This value coincides with that provided by 
the slow-roll formalism of inflation^ ^ ^ . 




3.2. Number of e-folds during inflation 

Another point we investigated refers to the number of e-folds during which inflation 
occurs, but for that we required the help of numerical solutions. 

First of all, we have to set initial conditions for 0j, 0j and Hi, which should be 
translated into the initial conditions Xi, yi and Zi. The different physical constraints 
to be imposed upon them are discussed in the following list. 



The first constraint to be taken into account is (j)i/mpi ^ 0{1). In order to 
obtain results as general as possible, we took the inequality 

<10. (17) 




mpi 

• The second constraint arises from the requirement that p^^i < rripi , a con- 
straint that is called the quantum boundary in Ref.^. With the help of the 
Friedmann constraint ([3]) , one finds that such constraint corresponds to the 
lower bound Zi > m/mpi, which is in turn equivalent to impose Hi < mpi. 

• The initial value (pi is found in terms of the Friedmann constraint xf = 
1 — yf, once yt is found from the above constraints. 

• Last constraint refers to inflationary solutions, and then Zi < Zend- 

We noticed, in the diverse numerical experiments we performed, that the number 
of inflationary e-folds mainly depends upon the initial value Zi. This is because the 
dynamical system reaches the attractor points very quickly, and then its subsequent 
evolution can be well described by the motion of the attractor points themselves. 

Thus, a formula for the number of inflationary e-folds is given by a combination 
of Eqs. ([Sc]) and (fT3|) . namely, 

N{z-Zi) = \n(-]^- X , (18) 
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Needless to say, we have found very good agreement between our numerical results 
and the analytic ones provided by Eq. (|18p . 

The total number of e-folds in a given inflationary stage, defined through 
Ntotai = N{zend',Zi) in Eq. (jlSp. is uniquely determined by the initial value Zi, 
or equivalently, by the initial value of the Hubble parameter Hi. If the total num- 
ber of e-folds should at least be Ntotai > 60, then the initial condition should be 
Zi < 0.15; this guarantees that there is sufficient inflation whatever the initial con- 
ditions on the scalar field are, as long as they are in agreement with the Friedmann 
constraint. 

The attractor solution indicates that the end of the accelerating stage happens 
once the Hubble parameter is Hend = m/-\/2, and the oscillations begin a bit later, 
when Hose = 2m/3, once the attractor point disappears. These two events are 
almost instantaneous, as the number of e-folds in between, as calculated again from 
Eq. nil), is AN = 0.05. 

3.3. Inflationary quantities 

We give here some instances about how we can write different inflationary quantities 
of physical interest in terms of the attractor values of our dynamical variables. 

The flrst one is the amplitude of primordial quantum perturbations of the infla- 
ton field given by 



where it is implicitly assumed, as usual, that all quantities on the r.h.s. are evaluated 
at the time the corresponding wave number k leaves the horizon, i.e. k = aH. Also, 
notice that the last term on the r.h.s. should be calculated using the attractor 
solutions in Eq. p^ . 

The last statement is not trivial, as the standard procedure is to calculate the 
values of H and (p in the slow-roll regime. This is not our case, as we are giving 
the aforementioned dynamical variables their values at the attractor points, and the 
latter are not subjected to any order approximation. The same remark applies to 
other quantities shown below. 

The wave number k in Eq. (jl9p can be given in terms of the smallest scale 
that leaves the horizon at the end of inflation, which we call kend- In terms of our 
variables, it is easy to show that 



That given, we can calculate the spectral index n{k) with the help of Eqs. (fig]) 
and (|20p . and then we obtain 




(19) 



ln(fc/fce„d) = N{z; Zi) - Ntotai + ln{Zend/z) . 



(20) 



n{k) - 1 



dlndl 



(k) 



2[z2-3x2(z)] 



(21) 



din fc 



[l-3a;2(z)] [l-2:r2(z)]' 
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where, as before, the last term on the r.h.s. is the final expression found from the 
attractor solution in terms of our dynamical variables. 

As said before, our calculations rely on exact inflationary attractor solutions, 
so that we can proceed further and calculate higher order derivatives of the scalar 
power spectrum. For instance, the running of the spectral index a can be obtained 
just by one more parametric derivation of Eq. |[2T|) . 

dn, ^ Axl{z)z^ [z2 + 9x^(z)-6] 
"~dlnfc 3[l-3a;2(z)]^[l-2a;2(z)]3 ■ 

As in other results shown before, this formula is evaluated at the attractor solution. 
By contrast, the slow- roll formalism needs of a second order calculation in the slow- 
roll parameters to get a precise result^ l-Sl, 

For completeness, we give below the values of the aforementioned inflat ion ary 
quantities at a time corresponding to 60 e-folds before the end of inflatiorff^. In 
terms of our control parameter, that time corresponds to z ~ 0.15, as ^^(-^2; 0.15) = 
60, see Eq. p^ . The resulting values are 

(51,''°^ -3.0105x102^, (23a) 



m 



2 

PI 



~ 0.969659 , (23b) 
(60) ^ „4.65 X lO^-* . (23c) 



The above values can be compared to those inferred from observations of the Cos- 
mic Microwave Background. For example, from the measured amplitude 5h — 
1.9 X 10^^ one obta ins the known result for the scalar field mass, which is 



m 



4. Conclusions 

We have shown, with the help of the theory of dynamical systems, the existence of 
attractor solutions in a scalar field theory minimally coupled to gravity and endowed 
with a quadratic potential, the simplest and most known of the monomial potentials 
in inflationary theory. 

The attractor solutions and other subsequent results were given explicitly in 
terms of the auxiliar variable z = m/H, which acted as a free parameter and con- 
trolled the behavior of the critical points. We called this method the fixed point (FP) 
approximation. In particular, we applied the formalism to the study of inflationary 
dynamics with a quadratic potential. Our results are explicitly based on attractor 
solutions, and then we were able to calculate diverse inflationary quantities without 
the need of order approximations. 

Also, we would like to stress out that the conditions for an inflationary stage 
are largely determined by the value of our control variable z, as it alone determines 



^See Ref.U^and table of parameters in |http:/ /lambda. gscf.nasa.gov| 
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the existence of an inflationary attractor solution. This confirms the results of pre- 
vious studies about the feasibility of inflationary solutions from arbitrary initial 
conditionJ^, see also the discussion in Refs. [20l [21] 

We believe that the FP approximation presented in this paper can be used 
further to fully exploit the predictions of the simplest chaotic inflatio nary model, 
and to have a more precise comparison with cosmological observationJ^SESIllllSl 
The FP approximation can be easily generalized to other inflationary model^; this 
is work in progress that will be published elsewhere. 
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Appendix A. Comparison with slow-roll 

In this appendix, we are going to make the comparison of the 'fixed point' (FP) 
aproximation to the slow-roll (SR) one. The SR approximation for the equations of 
r rLotion of a scalar field in a FRW metric and endowed with a quadratic potential 
jjnUni ^gee also Refs. [21 [Sjfor an earlier study of the SR equations) 

^^_!^, H^^^-^m^c^\ (A.l) 
^ 3H ' 6 ^ ^ ' 

In terms of our auxiliary variables (O, equivalent expressions for the SR equations 

are 

xsR--^, vIr-1- (A.2) 

From Eqs. (|13p . it can be seen that the SR solution is a first order approximation 
to the FP one. Actually, in the limit z <C 1 we get 

=^ t| , ys^±{l~'j^ ■ (A.3) 

The second equation shows that the correction to the usual SR assumption ysR — 1 
is of the second order. Likewise, the values of the inflationary quantities in the limit 
z < 1 are, from Eqs. (HH), JUj), and (|22|) . 

'-^W.^r^ n(fc)-l.--z^ a^--z^. (A.4) 



During the revision process of this work, we discovered that the FP approximation was success- 
fully applied to a quartic monomial potential in Ref. |24| where the authors also make a detailed 
comparison with the results obtained from slow-roll. 
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The above formulas can be written in terms of the number of e-folds given by 
Eq. (HH); to first order we find N{^/2, Zi) ~ {i/2)z^'^ in the hmit Zi < 1. If we use 
this last result in Eqs. (|A.3p . we find the usual expressions of the SR approxima tion 
for the different inflationary quantities in terms of the number of e-folds^j; see alsc|24l. 

Some remarks are in turn. First one is that the Friedmann constraint is not fully 
accomplished in the SR approximation, as we have seen that the deviations are of 
second order. Second, it is manifest that SR equations (|A.2p are not solutions of the 
original equations of motion. Moreover, the SR points are not even critical points 
of the dynamical system ([5]), and so it is difficult to study their attractor properties 
using standard techniques of dynamical systems. 

The SR and the FP trajectories coincide in the very early universe, in the limit 
z <C 1. But even for moderate small values of z, the deviations of SR from FP start 
to be significant, see Fig. [21 and then one has to include correction terms in the SR 
approximation at different 

levelPini. 
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